Cavity Optomechanics with a Bose-Einstein Condensate: Normal Mode Splitting 
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We study the normal mode splitting in a system consisting of a Bose-Einstein condensates (BEC) 
trapped inside a Fabry-Perot cavity driven by single mode laser field. We analyze the variations in 
frequency and damping rate of collective density excitation of BEC imparted by optical field. We 
study the occurrence of normal mode splitting which appears as consequences of hybridization of 
the fluctuations of intracavity field and condensate mode. It is shown that normal mode splitting 
vanishes for weak coupling between condensate mode and intracavity field. Moreover, we investigate 
the normal mode splitting in the transmission spectrum of cavity field. 
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I. INTRODUCTION 

Nano optomechanical systems that couple optical de- 
gree of freedom to the mechanical motion of a cantilever 
have been extensively investigated in recent past [1]. In 
such systems, coupling is obtained via radiation pressure 
inside a cavity [2-4] , and indirectly via quantum dots [5] 
or ions [6]. Significant progress has been made in the 
investigation of optomechanics, such as squeezing [7, 8], 
ultrahigh precision displacement detection [9-11], mass 
detection [12], gravational wave detection [13, 14], and 
the transition between classical and quantum behaviour 
of a mechanical system [15]. Moreover, entangling the 
electromagnetic field with motional degree of freedom of 
mechanical systems have been explored in various ap- 
proaches [16, 17]. The entanglement generated in such 
systems is significant both philosophically as well as tech- 
nically in relation with quantum informatics [18]. Fur- 
ther advances in experimental techniques make possible 
to couple the mechanical resonators with the statistical 
ensemble of atoms. In such systems, the interaction is 
mediated by the field inside the cavity that couples the 
mechanical resonators to the internal and motional de- 
grees of freedom of the atoms [19-22]. 

The normal mode splitting (NMS) is the coupling of 
two degenerate modes with energy exchange taking place 
on a time scale faster then the decoherence rate. More- 
over, the NMS is a phenomena ubiquitous in both the 
classical and quantum physics. The NMS, in the trans- 
mission spectrum of the cavity filed, has been observed 
when atoms are coupled to cavity field [23]. Moreover, 
the NMS has also been studied with artificial atoms in 
circuit quantum electrodynamics (QED) [24] as well as 
in single quantum dot cavity QED [25]. In this paper, we 
study the occurrence of the NMS through a new optome- 
chanical system, i.e, a coupled BEC-cavity configuration 
[26-34] . Briefly, the system consists of Fabry-Perot cav- 
ity containing a condensate that effectively behaves like 
a vibrating mirror [28]. In such systems, the collective 
density excitation of the BEC serves the analogy of a 
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moving mirror coupled to cavity field via radiation pres- 
sure force. The strong coupling of the quantized cavity 
field with the collective oscillations of the BEC has been 
experimentally analyzed with all the atoms are being in 
the same motional state [27, 28]. Our study reveals that 
the frequency and the decoherence rate are sensitive to 
radiation pressure force. We further analyze the occur- 
rence of NMS in the position spectrum of the conden- 
sate mode. It is observed that the position spectrum of 
the condensate mode splits into two peaks when coupling 
between condensate mode and intracavity field is consid- 
ered. In addition, we also discuss the NMS in the fluctu- 
ations spectrum of the light field emitted by the cavity. 
We show that the distance between two peaks increases 
linearly with BEC-cavity field interaction. The paper is 
structured as follows: In Section 2, we give the theoret- 
ical model of the system and its coupling with the en- 
vironment by using the Heiscnbcrg-Langevin equations. 
In section 3, we solve the dynamics in frequency domain 
and discuss the results. Finally, we provide concluding 
remarks in section 4. 



II. MODEL AND HAMILTONIAN OF THE 
SYSTEM 

We consider an ensemble of N two level bosonic atoms 
with resonant frequency u>t r is trapped inside a Fabry- 
Perot cavity and interact with standing- wave light field as 
shown in Fig.l. We assume the atom-field detuning A a is 
very large, therefore, one can adiabatically eliminate the 
excited atomic level. In the rotating frame at the driving 
field frequency ujp , the Hamiltonian of the Bose-Einstein 
condensate in the limit of small atom-atom interaction 
and small value of the external potential can be written 
as follows [28]: 

H = hA c Jc+hu> a ata-\ — -y=hg (a 1 "+a)c 1 ' c-ihE{c-c^) (1) 
v2 

where uj a = Auj r , A c = uj c -uj p + NU a /2, U = g 2 /A a , 
A a = ujp — uJtr, and g = y/~NU /2. Here, N is the num- 
ber of BEC atoms, u>t r is the transition frequency and 
g Q accounts for the coupling strength between the sin- 
gle intracavity photon and single condensate atom. The 
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FIG. 1. (Color online) A sample of two-level bosonic atoms 
with resonant frequency uj tr trapped inside a Fabry-Perot cav- 
ity of length L is interacting with standing laser field. Here, 
the left-end mirror is transmissive while the right-end mirror 
is perfectly reflecting. Moreover, the cavity is being driven by 
a laser of frequency uj v . 



first term in the Eq.(l) stands for the energy of the intra- 
cavity mode with creation (annihilation) operator c(c^) 
and frequency lo c . Moreover, the empty cavity resonance 
frequency w c is shifted due to the presence of the BEC in- 
side the cavity by an amount of NU /2. The second term 
describes the energy of the Bogoliubov mode of the col- 
lective oscillations of the BEC. Furthermore, a(a^) and 
uj a are, respectively, the creation (annihilation) opera- 
tors and frequency of the condensate mode. The dy- 
namics of the BEC can be explained as follows: The 
zero momentum state is only coupled to the symmet- 
ric momentum states, ±2?ifc, due to the absorption and 
stimulated emission of the cavity photons [28]. This can 
be explained as the condensate mode oscillates at fre- 
quency ui a = 4u> r = hk 2 /2m, where m is the mass of the 
atom and w r is recoil frequency. The third term describes 
the interaction between intracavity field with condensate 
mode. Moreover, g accounts for the strength of inter- 
action between field and BEC and it is clear from the 
expression of the g that the single atom-photon coupling 
is increased by the square root of the number of conden- 
sate atoms. The last term corresponds to the coupling 
between cavity mode and input laser field with coupling 
strength E and it is related to the input power P with 
\E\ = ^/2nP/hu! p , where n is the decay rate of the cavity 
field. 



condensate mode and intracavity field are given by 

a = -(iu a + l)a - i—f= c)c+ y/2rfa in , 
V2 

c = —(K + iA c )c — i—=(a + a^)c-\-E 
V2 



2 KjCi' 



(2) 



We assume that the noise associated with the light field 
is uncorrelated with the noise accounts for the conden- 
sate mode. For laser field, the noise and damping are due 
to the vacuum noise, losses from the cavity and fluctua- 
tions of the laser field. Moreover, the noise and damping 
accounts for condensate mode are due to the condensate 
atoms with additional nearby non-condensed atoms. In 
addition, a n and a,i n are the non-commuting noise oper- 
ators associated with optical field and condensate mode 
respectively. They have zero mean values and nonzero 
correlation functions [35]: 



(dal(t)da m (t')) 
(da m (t)dai(t>)) 
(dcl(t)dc m (t')) 
(dc ln (t)dcl(t')) 



n a 5(t-t'), 

(n a + l)S(t-t'), 

n c S(t-t'), 

(n c + l)S(t-t'), 



(3) 



where n c and n a are the occupation numbers of the op- 
tical and condensate modes, respectively. Moreover, all 
other correlations are zero. As we are in optical regime 
and a BEC at a temperature of at most a few //K, there- 
fore one can take n ca 0. 



III. DYNAMICS OF SMALL FLUCTUATIONS: 
NORMAL MODE SPLITTING 

In the following we linearized the operators in Eq.(2) 
around the steady state values, a = (a) , s +<9a, c = (c) ss + 
dc. Here, we have assumed that the fluctuation operators 
da and dc have zero mean. The steady state value of the 
intracavity mode is (c) = E/(k + i A), where the total 
effective detuning is 



In order to describe the complete dynamics of the 
subsystems involved in this problem, an adequate choice 
is to use the formalism of the quantum Langevin equa- 
tions. According to the Hcisenberg-Langcvin equation 
of motion, the commutation relations [a, a'] = 1 and 
[c, c^} = 1, the time evolution of the c and a, can be 
obtained. We derive the Heisenberg-Langevin equations 
for canonical variables and introduce the noise operators 
Ci n {t) and a,i„ weighted with the rates k and 7 which 
describe the dissipation of the intracavity field and con- 
densate mode of collective oscillations of the BEC re- 
spectively Therefore, Heisenberg-Langevin equations for 



For the sake of simplicity we assume that the field is real 
positive and this can be achieved by adjusting the phase 
of the laser field. Similarly, the steady state value of the 
condensate mode is (a) ss = [— ig ac /V2 (7 + iw a ] (c) ss - 
We linearize the Langevin equations of motion given in 
Eq.(2), and assume that pump field is intense and keep 
terms only up to first order in the fluctuation opera- 
tors. We rewrite each Heiscnberg operator in Eq.(2) as a 
sum of steady state value and fluctuation operator with 
zero mean value. Therefore, the linearized Heisenberg- 
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FIG. 2. (Color online) (a) Plot normalized effective frequency 
Weff/coa of the condensate mode as a function of the nor- 
malized frequency uj/ui a . Parameters values are ui a = 2n 4 x 
3.8 kHz, 7 = 27rx0.4kHz, A = u a , k = 0.1 uj a and G = 0.1 w a 
(solid blue curve), G = 0.2 ui a (dashed green curve), and 
G — 0.3 uj a (red dotted line), (b) Plot of normalized effective 
damping rate r y e ff/cj a of the condensate versus normalized 
frequency u)/u a for three different values of the BEC-field 
coupling, G — 0.1 uj a (solid blue line), G — 0.2 uj a (dashed 
green line), and G — 0.3 ui a (red dotted line). 



Langcvin equations are, 



da 



,G 



(iuj a + 7)9a - i— (dc + 9c 1 ") + y/5y a v 



,G 



dc = — (i A + «) dc — i— (da + <9a 1 ) + v / 2kc» t j(.5) 

The linearized quantum Langevin equations show the 
fluctuations of the Bogoliubov mode as the collective os- 
cillation of the BEC. The condensate mode is now cou- 
pled to the cavity field quadrature fluctuations by the 
effective couplings G = \[2gc s which can be made very 
large by increasing the amplitude c s of the intracavity 
field. Linearized quantum Langevin equations (5) and 
their corresponding Hermitian conjugate form a system 
of four first-order coupled operator equations, for which 
the Ruth-Hurwitz criteria [36] implies that the system 
will be stable only if the following stability conditions 
arc satisfied, 

51 = 2ft 7 {[ft 2 + K-A) 2 ][ft 2 + (w a + A) 2 ]+ 7 [(7 + 2ft) 

x(ft 2 + A 2 ) + 2kw 2 ]} + Acj q G 2 (7 + 2k) 2 > 

5 2 = Lj a (K 2 + A 2 ) -G 2 A > 0. (6) 

We choose a parameter regime such that the above 
stability conditions are satisfied to ensure stability of 
the system. The Hcisenbcrg-Langevin equations (5) 
are linear in creation and annihilation operator. We 
define the quadratures, i.e., dx a = (da + 9a t )/v / 2, 
dp a = (da — dat)/i*j2, dx c = (dc + <9ct)/\/2 and 
dp c = {dc — 9ct)/i\/2, and solve the Langevin equa- 
tions in Fourier space. Here, we introduce the dis- 
placement spectrum of Bogoliubov mode of collective os- 
cillation of the BEC which is obtained from the two- 
frequency auto-correlation function (dx a (ui)dx a (uj')) = 
S Xa (uj)S(uj +lu'). Therefore, the displacement spectrum 



S Xa (uj) = (1/2tt) / e-^-^'^dxa^dxa^^duj' of the 
condensate in Fourier space is given by 



S Xa (u) 



|d( W )|s 



2G 2 ft(ft 2 + Lo 2 a + A 2 ) 
(ft 2 + A 2 - w 2 ) 2 + 4k 2 w 2 



27 



, (7) 



where uj stands for frequency and 



d(ui) 



■AG 2 



«W7 ■ 



(ft 



is the modified susceptibility of the condensate mode due 
to radiation pressure. The effective frequency (w e //) of 
collective oscillation of the BEC and effective damping 
rate ("feff) of collective oscillation of the BEC are given 
by 



^eff = 
7e// = 



7 + 



w q AG 2 (k 2 — uj 2 



A 2 



(>< 2 



■jj 2 + A 2 ) 2 
2w a ftAG 2 



(2kw) 



(ft 2 



+ A 2 ) 2 + (2ftw) 5 



(8) 
(9) 



The frequency of the condensate mode is modified due to 
radiation pressure as shown in Eq. (8) and this is equiva- 
lent to the optical spring effect in case of opto-mcchanical 
system with moving mirror. The spectrum of conden- 
sate mode is described by the effective susceptibility d(oj) 
and displacement spectrum of collective oscillation of the 
BEC which consists of two terms, the first term is propor- 
tional to quantum fluctuation of the radiation pressure 
and, second term arises from quantum noises associated 
with the matter waves. Therefore, the position spectrum 
of the condensate mode is determined by radiation pres- 
sure and quantum noise. 

For our numerical calculations, we choose the pa- 
rameters very close to the BEC-cavity system [28, 31]. 
The recoil frequency ui r = 2ir x 3.8 KHz, decay rate 
7 = 2tt x 0.4 KHz, atom-field detuning A a = 2tt x 32 GHz 
and single atom-photon coupling is g = 2tt x 10.9 MHz. 
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FIG. 3. (Color online) The displacement spectrum of the 
Bogoliubov mode S Xa (w) as a function of the normalized de- 
tuning A/w a and normalized frequency uj/uj a for k = 0.1 uj a 
and (a) G — 0.1 u) a (b) G — 0.2 u a . The other parameters 
values are the same as in Fig. 2. 
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The optical spring effect leads to the shift in the fre- 
quency of the condensate mode and this shift is small 
for low-frequency oscillators [37] and shift is significant 
for high-frequency oscillators [38]. In Fig. 2(a), we plot 
normalized effective frequency u) e ff/u> a of Bogoliubov 
mode of collective oscillation of the BEC as a function 
of the normalized frequency u> /uj a for three different val- 
ues of the BEC-intracavity field interaction, G = 0.1 uj a 
(solid blue curve), G = 0.2 u a (dashed green curve), and 
G = 0.3 uj a (red dotted curve). It is noted that the devi- 
ation in the bare frequency of the condensate mode u> a is 
increased as the coupling between BEC and intracavity 
field is increased. In Fig. 2(b), we plot the normalized ef- 
fective damping rate 7 e ///w Q of the condensate mode ver- 
sus normalized frequency u> /to a for different values of the 
BEC-ficld interaction. The effective damping rate of the 
condensate mode is significantly increased as the BEC- 
field interaction is increased as shown in Fig. 2(b). There- 
fore, strong atom-field interaction causes higher atomic 
loss and this atomic loss is maximum near the resonance. 
The back action of the light causes the heating of the 
atoms, and as a consequence, atomic loss was observed 
in [39] . The increase in the effective damping is at the ba- 
sis of cooling of Bogoliubov mode of the BEC. The BEC 
couples to the intracavity field due to radiation pressure 
force and radiation pressure behaves as an thermal reser- 
voir for the BEC resonator. The effective temperature of 
the Bogoliubov mode of the condensate takes the value 
between initial thermal reservoir temperature and opti- 
cal reservoir. Therefore, when G > 7 one achieves the 
ground state cooling of the Bogoliubov mode. 

In Fig. 3, we plot the displacement spectrum of the 
BEC as a function of the normalized detuning A /u a and 
normalized frequency w/cj q for two different values of 
the BEC field interaction, i.e, G = 0.1 cj a , Fig. 3(a), and 
G = 0.2 oj a , Fig. 3(b). For small value of interaction be- 
tween condensate mode and cavity field, i.e G = 0.1 oj a , 
the fluctuations spectrum of the Bogoliubov mode barely 
shows the normal mode splitting. As we increase the 
coupling strength, the normal mode splitting becomes 
observable. For G = 0.2w a , the normal mode splitting is 
clearly seen in Fig. 3(b). It is observed that the normal 
mode splits in the presence of the large coupling between 
BEC and field. The normal mode splitting is due to the 
hybridization of the two oscillators: the electromagnetic 
field and collective oscillations of the BEC mode. The 
normal modes splitting results from the mixing of the 
fluctuations of the cavity field with the fluctuations of 
the condensate. The normal mode occurs at frequen- 
cies c4 = (A 2 + uj 2 a ± ^(A 2 -w 2 ) 2 +4G 2 w a A)/2. In 
this expression, we neglect the cavity decay rate re and 
damping rate of the Bogoliubov mode. The splitting is 
directly proportional to the coupling rate between Bo- 
goliubov mode and intracavity field, i.e, lj + — w_ oc G. 

The atom-field system is probed by measuring the 
optical field transmitted by the cavity, particularly, its 
spectrum. The power spectrum of the light emit- 
ted by the cavity is determined from Eq.(5). In fre- 
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FIG. 4. (Color online) The emission spectrum S C out(iS) as a 
function of the normalized frequency ui/uj a for different values 
of the BEC-field coupling strength, G = 0.1 u a (solid blue 
curve), G = 0.2u a (dotted red curve) and G = 0.3 u) a (green 
dashed curve). Moreover, the cavity detuning A = u a and 
the other parameters are the same as in Fig. (2). 



quency space, the fluctuating dc(uj) of the intracav- 
ity field can be obtained by using Heisenberg Langevin 
equations (5). With the quantum input-output rela- 
tion [40] dc out = \/2re<9c(aj) — dci n (uj), the power spec- 
tral density of the cavity output field is defined by tak- 
ing the Fourier transform of its autocorrelation func- 
tion {c\ ut {u}') c out (u)) = S cout {uj)S(uj' +w) as S cout (uj) = 
(1/2tt) / e- u ^^(dc cou t(uj)dc cou t(io'))duj. Therefore, 
the power spectra of the light emitted by the cavity is 
given by, 



S C out{^>) — 



[4re 2 |a(u;)| 2 + 2 7 |/3H| 2 ], (10) 



where, 



WaG 2 



a(u)) 
I3(lu) 

(1 c (uj) = (re - iuj) 2 + A 2 - 



V2(^ a - u 

— iV^KLUaG 



2 — iuj^f) 
K + iA 



V2AG 2 



It is clear from Eq.(10) that the spectrum of the output 
field consists of two terms: the first term describes the 
coupling of intracavity field with the Bogoliubov mode of 
collective density oscillations of the BEC and second term 
accounts for the vacuum noise associated with condensate 
mode. 

In Fig. 4, we plot the fluctuations spectrum of the 
output field as a function of the normalized frequency 
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u)/us a for different values of coupling strength between 
intracavity field with BEC. For small coupling between 
condensate mode and cavity field G = 0.1 oj a , one can see 
that only single peak (solid Blue curve) is appeared in 
the fluctuations spectrum of the output field. However, 
as the interaction of the intracavity field with BEC is 
G = 0.2 uj a increased, the spectrum splits into two peaks 
(Red dotted curve). For G = 0.3 w a , we observe S cou t(ui) 
further splits into two sideband peaks. This splitting is 
pure result of the coupling between condensate mode and 
cavity mode. It is also observed that the separation be- 
tween two peaks in the S cou t(uj) is directly proportional 
to the BEC-field interaction strength. This splitting de- 
pends linearly on the coupling parameter G and goes to 
zero in the absence of BEC-field coupling. 



IV. CONCLUSION 

In conclusion, we theoretically analyze the normal 
mode splitting in a system which consists of Bose- 
Einstein condensate trapped inside a Fabry-Perot cav- 
ity driven by laser field. We observe that the frequency 
and damping rate of the condensate mode display shift 
due to the radiation pressure force. The hybridization of 
the Bogoliubov mode with the fluctuations of the cavity 
field leads to normal mode splitting. The results show 
that normal mode splitting in fluctuations spectrum of 
the condensate mode depends linearly on the coupling 
of the BEC with optical field. We also analyze the cou- 
pled dynamics of the BEC and cavity field by probing 
the emission spectrum of the cavity. Normal mode split- 
ting is also observed in the transmission spectra by cavity 
and this splitting is directly proportional to the BEC-field 
coupling strength. Therefore, the coupling rate between 
Bogoliubov mode and interactivity field can be deter- 
mined by measuring the distance between two splitting 
peaks. 
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